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INTRODUCTION 


In  the  second  year  of  Contract  F 1 963-8 l-K-0012 ,  "Ground 
Response  Alluvial  Basins  Due  to  Seismic  Disturbances",  from 
March  26,  1982,  our  research  emphasis  has  been  In  the 

following  four  areas: 

I.  Solving  the  finite  element  basin  models  for  elastic 
(P  and  S)  waves,  models  considered  are  relatively  more  complex 
than  what  we  studied  In  the  period  of  the  first  year  of  the 
Contract . 

II.  Studying  the  finite  element  source  mechanism  for  the 
cases : 

(A)  One  Degree  of  Freedom  SH-wave  Problem  with  : 

(I)  concentrated  line  source, 

(II)  concentrated  coupled  line  source. 

(B)  Two  Degrees  of  Freedom  Elastic  Wave  Problems  with: 

(I)  directional  concentrated  line  source, 

(II)  concentrated  coupled  line  source. 

(III)  non-d 1 rec t 1 ona 1  line  source, 

III.  Modifying  and  Improving  our  two-dimensional  finite 
element  codes  for  the  cases  of  elastic  and  Sh-  waves  by: 

(A)  replacelng  the  4CST(Four  Constant  Strain  Triangles) 
algorithm  by  the  Averaging  2CST  algorithm, 

(B)  replacing  the  element-oriented  approach  by  the 
noda 1 -po 1 nt-or i ented  approach, 

(C)  changing  the  structure  of  the  global  assemblage 
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>  system. 
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I  IV.  Implementing  the  two-dimensional  finite  element 

codes  with  a  snapshots  plotting  software. 
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SCIENTIFIC  ACCOMPLISHMENTS 


I.  Finite  Element  Modeling  of  Preliminary  Basins 

(A)  Partial  Basins  with  Display  of  Synthetic  Seismogram 

As  a  continuation,  we  have  studied  two  preliminary  basin 
models  as  shown  In  Figures  1  and  2  for  the  case  of  elastic  (P 
and  S)  waves.  The  results  for  the  corresponding  SH-wave  case 
have  been  reported  In  the  Scientific  Report  No.  1  (Kuo  and 
Teng,  1982).  In  Figure  1  (Model  3-A  In  the  Scientific 
Report  No.  1),  a  homogeneous  material  with 

Vp  ■  5,000  m/sec ,  ■  2,900  m/sec,  ■  2.67  gm/cm3  ,  is 

excited  by  a  vertical  function  FQ(t)  as  shown  In  Figure  7-B 
at  Point  S  on  the  free  surface  of  the  deeper  basin  DC,  while 
the  receivers  are  located  on  the  free  surface  of  the  shallower 
basin  In  the  region  BC .  In  Figure  2  (Model  3-B  In  the 
Scientific  Repeort  No.  1,  Kuo  and  Teng,  1982),  a  vertical 
source  Is  located  In  the  shallower  basin,  while  the  receivers 
are  In  the  deeper  basin.  Again,  the  forcing  function  Is  the 
first  derivative  of  a  Gaussian  function  with  Its  center 
frequency  f c  ^  2.7  Hz.  Figures  3  and  4  are  the  results  of 
the  displacement  components  along  plane  ABC  In  both  models 
respect  1 ve ly . 


DISPLACEMENT 


(B)  Partial  Basins  with  Display  of  Snapshots 

In  order  to  obtain  clearer  display  of  the  results  for  the 
basin  modeling,  an  Aldridge  Laboratory  computer  code  for 
snapshots  has  been  Implemented  to  our  finite  elemenet  computer 
codes.  The  display  of  snapshots  provides  a  technique  whereby 
the  time  history  of  the  Interaction  can  be  visualized  and  the 
resulting  wave  shapes,  direction  of  propagation  and  nature  of 
mode  conversion  can  be  Identified.  For  demonstration,  the 
snapshots  have  been  made  for  the  following  three  cases: 

Case  I:  Model  3-A  for  Elastic  Waves 

Figures  5-A  to  5-G  are  snapshots  of  the  synthetic 
wavefleld  at  0.3  sec  Intervals.  Vertical  and  horizontal 
displacement  components  are  given  separately.  Each  snapshot 
Is  Independently  scaled  and  the  amplitudes  magnified  10  times. 
Black  Is  positive  displacement  and  white  Is  negative 
displacement.  Diffracted  waves  at  the  square  corners  are 
clearly  Identified  after  0.9  sec. 


Figure  5-G .  A  Snapshot  View  < 
after  the  Detonation  of  the 
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Case  2:  Model  3-A  for  SH-Wave 

In  Figure  6,  a  SH-wave  sourcj  F^  ( t )  Is  located  at  the 
corner  S  of  1  he  lower  basin.  The  scurce  function  F^(t)  Is 
shown  In  F1gur€  7-A.  Figures  8-A  to  8-0  are  the  snapshots  of 
the  synthetic  wavefleld  at  0.2  sec  Intervals.  In  this  simple 
one-degree-of-f reedom  problem,  the  reflections  from  the  free 
surfaces  and  the  diffractions  from  the  corners  can  be  clearly 
identified.  Figures  8-M  to  8-0,  we  observed  that  the 
waveflelds  are  contaminated  by  the  undesired  reflections  from 
the  side  and  the  bottom  boundaries. 


Figure  7-B.  Source  Function  for- Elastic  Cases. 


the  ion  >>f  th 


km 


Snapshot  View  of  Have  Pattern  at 
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Case  3:  Model  4  for  SH-Wave 

Figure  9  Is  the  geometry  of  a  preliminary  basin  model 
with  an  Inclined  slope  of  45  .  A  SH-wave  source  function 

( t )  (Figure  7-A)  is  located  at  the  corner  S  .  Figures 
10-A  to  10-0  are  the  snapshots  of  the  synthetic  wavefleld  at 
0.2  sec  Intervals.  Two  seconds  after  detonation  of  the  source 
undesired  ref  1  actions  start.  Comparing  the  results  of  this 
case  with  these  of  Case  2,  a  different  phenomena  can  be  seen 
In  Figures  10-'<  to  10-0.  In  the  neighborhood  of  the  source 
location,  the  medium  is  quiescent  2.2  sec  after  detonation  of 
the  source.  In  Case  2,  however,  there  Is  no  such  quiescent 
period.  A  reflection  appears  from  the  vertical  portion  of  the 
free  surface  boundaries. 


I 


the  Detonation  of  the  Source. 
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II.  Finite  element  Source  Mechan1?m 


In  finite  element  algorithms,  the  simulation  of  a 
non-d 1 rect 1 ona 1  line  or  point  source,  (the  cylindrical  source 
for  two-dimensional  problem,  or  the  spherical  source  for 
three-dimensional  problem),  for  the  two-  or 
three-deg ree-of -f reedom  elastic  wave  case  has ;  been  a  very 
challenging  subject.  The  simplest  or  the  most  natural  way  to 
excite  the  elastic  medium  Is  to  apply  a  directional  forcing 
function  to  the  structure.  Likewise,  for  the 
one-degrea-of-f reedom  problem,  special  effort  Is  needed  to 
simulate  a  couple  source.  In  the  Scientific  Report  No.  1  (Kuo 
and  Teng,  1982)  of  the  present  Contract,  we  proposed  to 
Investigate  systematically  this  subject  by  using  a 
energy-sha r 1 ng-noda 1 -po 1 nts  technique.  During  the  past  year, 
we  have  studied  the  two-dimensional  source  mechanism  problems 
for  the  cases  of  SH-waves  and  elastic  waves  (both  P  and  S) 
respectively.  For  the  SH-wave  case,  two  types  of  sources  are 
considered:  (1)  concentrated  line  source,  and  (11) 
concentrated  coupled  line  source.  For  the  elastic  wave  case, 
three  types  of  sources  are  considered:  (1)  concentrated 
vertical  (directional)  line  source,  (11)  concentrated  coupled 
line  source,  and  (ill)  concentrated  non-d 1 rect 1 ona 1  line 
source . 


Consider  a  two-dimensional  whole  space 
with  vn  ■  5,000  m/sec,  vc  ■  2,900  m /sec, 

P  * 

Figure  11  shows  the  finite  element  mesh.  In 


elastic  medium, 
y  *2.67  gm/cm^  , 
each  dimension. 


30  elements  are  used.  Point  S(0,0)  is  the  location  of  the 


external  source.  Using  this  model,  we  studied  the  following 


cases : 


(A>  One  Degree  of  Freedom  SH-Wave  Problem 

(1)  Concentrated  Line  Source 

As  we  mentioned  previously,  this  Is  the  most  natural  way 
to  apply  an  external  source  without  making  any  effort.  At 
Point  S(0,J0f>  In  Figure  11,  a  forcing  function  F-^(t)  (Figure 
7-A)  Is  applied.  Figures  12-A  to  12-H  are  the  displacement 
field  at  0.1  sec  Intervals.  The  cy 1 1 ndr 1 ca 1 ly  spreading  wave 
propagation  can  be  observed  clearly. 


(11)  Concentrated  Coupled  Line  Source 

In  order  to  simulate  a  coupled  line  source  at  point 

S(0,0)  of  Figure  11,  we  set  the  points  S^(0,*0),  and 

S2<0,~0)  occupying  the  same  location  as  point  S  and  apply  the 

forcing  function  F  (t)  to  each  Individual  point  but  In 

1 

opposite  directions  as  shown  In  Figure  13.  In  other  words, 
the  displacement  fields  are  w^(0,+0,t)  and  W2(0,~0,t)  at  S^ 
and  S2  respectively.  At  the  source  points  In  Figure  13,  dot 
means  outward  from  the  plane  of  the  paper,  cross  means  Inward 
towards  the  plane  of  the  paper.  Points  A,  B,  C,  and  D  are  the 
nodal  points  surrounding  point  S  In  Figure  11.  Figure  14-A  to 
14-H  are  the  snapshots  of  the  synthetic  displacement  field 
w(x,y,t)  at  0.1  sec  Intervals  showing  the  radiation  pattern 
of  a  coupled  source. 


Figure  lU .  Snapshot  View  of  Wave  Pattern  at  0.1  sec  Interval  After 
The  Detonation  of  the  Concentrated  Couple  Source. 


(B)  Two  Degree  of  Freedom  Elastic  Wave  Problem 

(1)  Directional  Concentrated  Line  Source 

For  the  elastic  wave  case,  the  most  natural  way  of  a 
source  Is  the  directional  source.  In  Figure  15,  a  vertical 
line  source  F0(t)  as  shown  In  Figure  7-B  Is  located  at  point 
S.  Figures  16-A  to  16-1  are  the  snapshots  of  the  absolute 
values  of  the  displacement  resultants.  The  time  Interval 
between  snapshots  Is  #.0575  sec. 


Detonation  of  the  Directional  Line 


((1)  Concentrated  Coupled  Line  Source 

Figure  17  shows  the  mechanism  of  the  coupled  source  for 

the  elastic  wave  case.  Points  S  ,<0,+0)  and  S  (0,-fif) 

1  2 

occupy  the  same  location  and  are  subjected  to  a  moment 
forcing  function  Fq  ( t ) ,  with  the  same  magnitude  but  opposite 
directions.  The  radiation  patterns  of  the  displacement 
magnitude  are  given  In  Figures  18-A  to  18-H. 
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(111)  Non-D 1 rect lonal  Concentrated  Line  Source 

In  order  to  simulate  a  non-d 1 rect Iona  1  concentrated  line 
source  at  point  S(0,0)  In  Figure  11,  we  set  the  points 
S^-0,+0)  ,S2<+0,+0>  .S^+0,-0)  ,  and  S^-0,-0)  at  the 
same  location  as  point  S  and  apply  a  forcing  function  F^(t) 
to  each  Individual  point  In  a  manner  shown  in  Figure  19. 
The  cy 1 1 ndr leal ly  spreading  radiation  patterns  of  the 
displacement  magnitude  can  be  observed  In  Figure  20-A  to  20-J . 


In  the  preceding  five  cases,  the  non-reflecting 
boundaries  technique  are  used  on  all  four  boundaries. 
Without  using  the  non-ref lect fng  bourdarles 

undesired  reflections  would  start  at  the  15th  time  step,  0.3 
seconds  for  tho  SH-wave  cases  and  0.1725  seconds  for  the 
elastic  wave  cases,  after  the  detonation  of  the  source. 


III.  Modification  of  the  Finite  Element  Codes 

In  the  old  version  of  our  two-dimensional  finite  element 

codes,  (SH-wave  and  Elastic  Wave  Cases),  we  used  the  4-CST 

(Constant  Strain  Triangles)  quadrilateral  elements  In  order  to 

avoid  the  space-grid  skewness.  The  degree  of  freedom  of  the 

Internal  node,  such  as  node  5  In  Figure  21-A,  Is  statically 

condensed.  The  degrees  of  freedom  of  such  nodes  do  not  appear 

$ 

In  the  global  assemblage  equations  so  that  the  storage  during 
the  computat  on  Is  reduced  cons  1 der ab 1 y .  The  4-CST 
formulation  has  been  a  very  good  approach  to  obtain  solutions 
for  static  cases.  However,  In  the  algorithm  of  4-CST  for 
dynamic  case.  It  was  assumed  that  the  Internal  condensed  nodal 
points  of  each  quadr 1 1  a tera 1  element  were  subjected  to  no 
inertial  forces.  In  the  present  modified  version  of  these 
codes,  an  average  2-CST  algorithm  Is  adopted.  In  Figure  21-B, 
a  quadr 1 1 atera 1  Is  composed  of  two  triangles.  Two  distinct 
ways  of  subdividing  the  qu3dr 1 1 atera 1  are  given.  Therefore, 
the  assembled  quadrilateral  stiffness  matrix  can  be  obtained 
by  linear  Interpolation  without  any  Internal  condensation. 
Based  on  the  averaging  2-CST  formulation,  a  better 
approximation  of  the  solutions  can  be  obtained. 


(B)  Nodal -Point-Oriented  Approach 


In  the  old  version  of  our  two-dimensional  finite  element 
codes,  the  global  stiffness  matrix  Is  obtained  by  the 
assemblage  of  the  stiffness  matrix  of  each  Individual 
quadralateral  element  for  each  time  step  In  the  time 
Integration . 

In  the  present  modified  version  of  the  codes,  the  global 
matrix  Is  assembled,  only  the  non-iero  stiffness  matrix 
members  are  Involved,  for  each  nodal  point  before  the 
computlon  of  the  time  Integration. 

The  advantages  of  the  old  codes  are: 

(1)  The  Incore  storage  required  Is  very  little,  only 
the  products  of  the  global  stiffness  matrix  and 
displacements  are  stored. 

(11)  Less  mistakes  can  be  made  particularly  for  the 
source  mechanism  problems. 

The  advantages  of  the  new  codes  are: 

(1)  The  computing  time  Is  drastically  reduced 

since  local  assemblages  for  each  time  step  are 
avoided.  The  amount  of  the  reduced  time  depends 
upon  the  size,  the  property  of  each  individual 
problem.  For  a  model  of  4 0x40  elements,  the 
computing  time  using  the  new  version  program  Is 
only  1/5  of  that  using  the  old  version. 

(11)  The  disk  storage  for  problems  with  Irregular  element 
meshes  is  .  the  same  as  that  of  the  problem 
with  regular  size  element  mesh. 


; 

I 


In  both  versions, 
excitation  method  are 


the  2-CST  formulation  and  the  effective 
used.  The  restart  back-up  option  Is 
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also  valid  for  the  new  version  codes 


FUTURE  RESEARCH 


The  seismic  velocity  models  for  western  alluvial  basins 
have  been  studied  extensively  by  Battls  (1981).  He  derived  a 
first-order  model  for  a  “typical"  MX  deployment  alluvial  basin 
with  the  features: 

(1)  The  model  basin  Is  approximately  70  km  In  length 
with  a  width  of  17  km, 

(11)  Th  maximum  depth  to  bedrock  within  the  valley  Is 
2000  m. 

(Ill)  The  basin  Is  enclosed  by  mountain  ranges  which 
parallel  the  major  axis.  These  ranges  are  about 
15  km  wide  with  an  average  elevation  of  750  m  above 
the  base  floor. 

( 1 v  )  The  character  1 st 1c  seismic  velocity  profile  of  the 
typical  basin  Is  as  follows: 


Layer 

Thickness 

(meters) 

Depth 

(meters) 

(kmfsec) 

V. 

(km/sec) 

1 

10 

0 

0.9 

0.3 

2 

115 

10 

1.6 

0.  55 

3 

75 

125 

2.0 

0.8 

4 

200 

200 

2.25 

1.  1 

5 

400 

400 

6 

1.5 

6 

450 

800 

3.  0 

1.75 

7 

(to  basement) 

1250 

3.  5 

2.  0 

8 

“  *  * 
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(v)  The  typical  basin  In  the  deployment  region  Is 
symmetric. 


Ha  further  pointed  out  that  the  concept  of  a  "typical*  basin 
model  Is  misleading.  Each  basin  has  Its  own  significant 
variations  from  the  typical  one.  In  studying  the  seismic 
responses  In  alluvial  basins,  both  the  typical  and  variations 
on  models  must  be  considered. 

Therefore,  even  a  typical  valley  structure  Is  a  complex 
Intermeshing  of  materials  with  high  velocity  gradients  at 
least  In  two  and  probably  all  three  dimensions. 

In  the  past  two  years,  our  research  emphasis  of  the 
present  contract  has  been  on  the  two-dimensional  studies 
Including  preliminary  basin  modeling,  source  mechanism,  and 
computer  codes  Improvements  In  finite  elements.  In  the  coming 
year  of  the  contract,  the  research,  within  the  limit  of  the 
Aldridge  PRIME  75fiT  computer,  would  be: 

(A)  Extending  the  above  two-dimensional  studies  to  three 
dimensions. 

(B)  Simulating  the  two-dimensional  alluvial  basin  models 
with  as  closely  as  possible  to  the  first-order 
typical  basin  model  derived  by  Battls  or  the  models 
provided  by  the  Air  Force  Geophysics  Laboratory  by 
finite  element  method  with  different  type  of  sources. 
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